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ABSTRACT 

We report the discovery of the sixth known eclipsing double white dwarf (WD) system, 
SDSS J1152+0248, with a 2.3968 + 0.0003 h orbital period, in data from the Kepler Mis¬ 
sion’s K2 continuation. Analysing and modelling the K2 data together with ground-based fast 
photometry, spectroscopy, and radial-velocity measurements, we determine that the primary 
is aDA-type WD with mass M\ = 0.47 + 0.11M©, radius R\ = 0.0197 + 0.00357?©, and cool¬ 
ing age t\ — 52 + 36Myr. No lines are detected, to within our sensitivity, from the secondary 
WD, but it is likely also of type DA. Its central surface brightness, as measured from the 
secondary eclipse, is 0.31 of the primary’s surface brightness. Its mass, radius, and cooling 
age, respectively, are M 2 = 0.44 + 0.09M©, 7?2 = 0.0223^ ^^7?©, and t 2 = 230+ 100 Myr. 
SDSS J1152+0248 is a near twin of the double-lined eclipsing WD system CSS 41177. 

Key words: methods: statistical - techniques: radial velocities - binaries: eclipsing - stars: 
individual: SDSS J115219.99+024814.4 - white dwarfs 


1 INTRODUCTION 

White dwarfs (WDs) in general, and WDs in binary systems in par¬ 
ticular, are important in a broad range of astrophysical contexts. To 
name some: for understanding the final stages of stellar and close¬ 
binary evolution (e.g. Benvenuto & Althaus 1998; Holberg et al. 
2012), as progenitors of hot subdwarfs (e.g. Han et al. 2002, 2003) 
and Type-la supernovae (e.g. Maoz et al. 2014); as sources of grav¬ 
itational waves (e.g. Nelemans 2009; Amaro-Seoane et al. 2013); 
and even as possible hosts of extrasolar planets on which life could 
potentially be detected (Agol 2011; Loeb & Maoz 2013). Eclipsing 
WD systems can be particularly revealing in that they permit direct 
measurements of the radii and masses of their constituents. 

Over the last years, the number of known binary systems with 
a WD component has increased significantly. There are more than 
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70 known double WDs (Marsh 2011; Kaplan et al. 2012; Brown 
et al. 2013; Gianninas et al. 2015). Since the 2010 discovery of 
the first eclipsing double WD binary, NLTT 11748 (Steinfadt et al. 
2010), four more such systems have been discovered (Parsons et al. 
2011; Vennes et al. 2011; Brown et al. 2011; Kilic et al. 2014). 
Among them, there is one known case of a double-lined eclipsing 
double WD system, CSS 41177, in which full model-independent 
derivation of the system parameters is possible (Bours et al. 2014). 

In this paper, we present the discovery of the sixth known 
eclipsing double WD system, SDSS J115219.99+024814.4 (here¬ 
after SDSS J1152+0248). SDSS J1152+0248 was identified as 
a g= 18.35 mag WD candidate in the Sloan Digital Sky Survey 
(SDSS) based on colour selection (Girven et al. 2011). It was in¬ 
cluded as a target for the 2-Wheeled Kepler continuation mission 
(K2) as a part of a program targeting -150 WDs (Kilic et al. 2013). 
As detailed below, photometry of the target, obtained during ‘Cam¬ 
paign 1’ of the K2 mission, revealed periodic (P = 2.4 h) primary 
and secondary eclipses. Follow-up spectroscopy confirmed the pri- 
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mary as a DA WD, and showed radial-velocity modulation with 
the same periodicity, while fast photometry revealed the details of 
the primary and secondary eclipses. Below, we describe these ob¬ 
servations, their analysis, and our modelling of the system and its 
parameters. 

2 OBSERVATIONS 

2.1 Photometry: K2 

SDSS J1152+0248 was observed in K2 Campaign 1, and was given 
the Ecliptic Plane Input Catalogue (EPIC) number 201649211. The 
campaign lasted from 2014 May 30 to 2014 August 21, and cov¬ 
ered a field in the North Galactic Cap. Kepler data are divided into 
‘cadences’. Our target’s data have a ‘long cadence’ in which 270 
frames are coadded on board, before the photometry of each target 
is downlinked to Earth, resulting in an integrated exposure time per 
epoch of 1766 s, or about half an hour. 1 

Analysis of the light curve, described in detail in Section 3.1, 
below, revealed periodic (2.4 h) primary and secondary eclipses, 
each of duration of 30 min in the phase-folded light curve. The pri¬ 
mary eclipse depth was ~ 2%. The similarity between the eclipse 
duration and the observing cadence suggested much deeper and 
briefer true eclipses, indicating a WD companion to the candidate 
primary WD. 

2.2 Spectroscopy: WHT 

A spectrum of the primary WD was obtained on 2015 March 16 
with the Intermediate-dispersion Spectrograph and Imaging Sys¬ 
tem (ISIS) on the 4.2-m William Herschel Telescope (WHT) in 
La Palma, Spain. A 1300 s exposure using the two channels of 
the instrument covered the ranges of 3600 to 5200 A and 5600 to 
8000 A at 0.8 A and 0.9 A resolution, respectively. This spectrum 
confirmed SDSS J1152+0248 as a DA-type WD (i.e. with a hydro¬ 
gen atmosphere). 

2.3 Photometry: Wise 

To investigate the suspected short eclipse duration, we obtained un¬ 
filtered imaging photometry using the PI camera on the 1-m tele¬ 
scope of the Wise Observatory in Israel, on the nights of the 2015 
March 17 and 22, covering about two full cycles of the system on 
each night. The observational cadences were 300 s and 30 s on the 
first and second night, respectively. These observations confirmed 
that the eclipse durations are of order 1 min, with primary and sec¬ 
ondary eclipse depths of about 50% and 10%, respectively. 

2.4 Spectroscopy: MMT 

Additional spectra of the system were obtained using the Blue 
Channel Spectrograph on the Multiple Mirror Telescope (MMT) 
in Mt. Hopkins, Arizona, to search for radial-velocity (RV) varia¬ 
tions in the primary WD, and for spectral signatures of the putative 
secondary WD. Thirteen exposures of 360 s each were taken with 
the 832 line mm -1 grating and a 1 arc sec slit on the night of 2015 
March 24, covering the range 3600 to 4500 A with 1.0 A resolution. 
As detailed in section 3.4, the spectra revealed RV variations at the 

1 See http://archive.stsci.edu/kepler/manuals/Data_Characteristics.pdf. 


expected period and phase, with semi-amplitude K\ « llOkms -1 . 
These spectra did not reveal any sign of the secondary, neither in 
the high-order Baimer lines (as might be expected if the secondary 
were another, bright enough, DA WD), nor in the presence of any 
other lines in the spectrum, to an equivalent-width limit of ~ 0.4 A. 
As a further test, we obtained additional spectra of the Her region 
of the spectrum, as the narrow non-LTE (NLTE) core present in 
the Her line of WDs permits resolving close spectral components. 
Four exposures, of 1080 s each, were taken with 1200 line mm -1 
grating and a 1 arcsec slit, on the night of 2015 March 27, covering 
the range 5855 to 7165 A with 1.5 A resolution, and were obtained 
approximately at the times of quadratures of the cycle (as predicted 
from our Kepler ephemeris for this system and from the RV curve 
obtained 3 nights before), in order to give the maximal spectral 
separation between the two components (if present). The spectrum 
reveals only a single Her component. However, based on the sys¬ 
tem parameters we derive below (see Section 3.5), we would not 
have expected to detect the Baimer lines from the secondary, given 
the signal-to-noise ratio of the spectrum. Although, in principle, the 
secondary could be a DC-type WD (a cool helium-atmosphere WD 
with a featureless optical continuum), this would be inconsistent 
with the lack of dilution we find in the primary’s Baimer lines (see 
Section 3.3). Furthermore, we will show that the secondary has a 
mass ~ 0.45M©, and Bergeron et al. (2011) have shown that helium 
atmosphere (DB-type) WDs rarely have masses < 0.5M©. Finally, 
we show in Section 3.5 that the results of our radial velocity esti¬ 
mates are likely affected by the presence in the spectrum of lines 
from the secondary WD. We conclude that the secondary is likely 
a DA WD, one whose lines and kinematics could be measured with 
deeper observations, which would then make this system the sec¬ 
ond known double-lined eclipsing WD system. 


2.5 Photometry: APO 

To probe the structure of the eclipses for the purpose of con¬ 
straining the system’s geometry, fast photometry was performed 
on the night of 2015 April 15 with the Agile photometer on the 
ARC 3.5-m telescope at Apache Point Observatory (APO) in New 
Mexico. A continuous sequence of 10 s exposures spanned two 
full cycles of the system. Observations used a BG40 filter, which 
is a broad-band filter covering ~ 3400- 6000 A. After dark sub¬ 
traction and flat correction, we retrieved aperture photometry on 
SDSS J1152+0248 and two nearby comparison stars. We chose 
aperture radii of 2 arcseconds which minimized the variance of the 
out-of-eclipse light curve. We constructed the light curve by divid¬ 
ing the flux from the target star by a linear combination of the fluxes 
of the comparison stars, and decorrelated trends with airmass. 


2.6 Photometry: McDonald 

Additional fast photometry was obtained on the night of 2015 April 
18 with the ProEM Camera on the 2.1-m Otto Struve telescope at 
the McDonald Observatory in Texas. Again, 10 s exposures with 
the BG40 filter (similar to the BG40 filter used at APO, see above) 
covered two full system cycles. We used the IRAF script ccd_hsp 
(Kanaan et al. 2002) to perform aperture photometry, and WQED 
(Thompson & Mullally 2013) to divide the target flux by the flux 
of a comparison star in order to correct for transparency variations, 
and to compute the B ary centric correction. 
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of its sample points. The error in each bin was estimated as 1.48 
times the median absolute deviation around the median, divided by 
y^bin, where ^bm is the number of points in that bin. 

Fig. 1 shows, for SDSS J1152+0248: the detrended K2 light 
curve; the Fourier amplitude spectrum, with a clear peak at 10 d -1 , 
and a strong harmonic at 20d _1 ; and the phase-folded light 
curve, revealing a primary and secondary eclipse, with a period of 
2.3968 + 0.0003 h. The second harmonic arises partly from the sec¬ 
ondary eclipse, which together with the primary eclipse introduces 
into the Fourier spectrum a double frequency, and partly from the 
highly non-sinusoidal form of the variations. The additional peaks 
in the Fourier spectrum result from higher harmonics, beyond the 
Nyquist frequency, that get reflected and shifted back into the spec¬ 
trum. As already noted, the eclipse duration of 30 min in the K2 
data, similar to the K2 time resolution per individual measurement, 
indicates a much briefer and deeper true eclipse that has been di¬ 
luted by the K2 cadence. 


Figure 1 . K2 data for SDSS J1152+0248. Top: reduced and detrended light 
curve; middle: FFT amplitude spectrum; bottom: phase-folded and binned 
light curve. Dashed red curve is best-fitting eclipse plus relativistic beaming 
model (see Section 3.2). 


3 DATA ANALYSIS AND MODELING 
3.1 K2 light curve and orbital period 

The periodic eclipses in SDSS J1152+0248 were discovered in the 
process of searching for periodicities in the K2 data for WDs. With¬ 
out two of its original reaction wheels, the images of the K2 mission 
are less stable than those of the original Kepler mission. 

Vanderburg & Johnson (2014) developed a self-flat-fielding 
photometry technique to reduce the effects of spacecraft roll on the 
K2 data. We used manually defined target apertures and the kepsff 
package under PyKE (Still & Barclay 2012) to reduce the K2 data. 
Self flat fielding reduced most of the spacecraft roll effects, but 
with remaining long-timescale variations. Additional artefacts in 
the data include small gaps of order an hour caused by the removal 
of low-quality cadences induced by maintenance procedures (such 
as angular momentum desaturation of the reaction wheels); a sin¬ 
gle ~ 3 day-long break in the data in the middle of the campaign, 
for downlinking collected science data to Earth; and an unevenly 
sampled grid that results from the conversion of Universal Time 
(UTC) to barycentre-corrected Julian date (BJD) in the middle of 
each cadence. 

To deal with these problems, we further processed each WD 
light curve through the following steps. The reduced light curve 
was reassigned to an evenly spaced grid using a nearest-neighbour 
interpolation, and gaps in the data were replaced by zeros. Post¬ 
reduction trends that remained in the light curve were removed by 
breaking up the data into 5 d intervals, and subtracting a linear fit 
from each interval, resulting in a trend-free light curve with zero 
mean. Data points deviating by more than 3 standard deviations 
from the mean were replaced by zeros. A Fourier amplitude spec¬ 
trum of the light curve was then calculated and used to get an ini¬ 
tial estimate of the period. The period was fine-tuned by inspecting 
the phase-folded light curve and adjusting the period until the least 
scatter was achieved. The scatter was determined by measuring the 
resulting modulation amplitude over a range of values close to the 
found period. The period errors were defined as the +cr intervals of 
the mean of a Gaussian fit of the scatter. Finally, the period-folded 
light curve was binned into 250 bins, each containing the median 


3.2 Eclipse modelling from fast photometry data 


We next combine the APO and McDonald fast photometry data 
to obtain time-resolved coverage of the primary and secondary 
eclipses. We then model these data to constrain some of the sys¬ 
tem’s geometric parameters. 

To obtain reliable estimates of photometric errors, each light 
curve, from APO and from McDonald, was assigned a constant er¬ 
ror of 1.48 times the median absolute deviation around the median 
of points outside of the eclipses. This resulted in errors of 5% for 
the APO data and 3% for McDonald. We tested for correlated errors 
among consecutive data points by calculating the auto-correlation 
function (ACF) for each of the light curves. Both ACFs fall to near 
zero in one time step, showing that any correlated errors are at a 
level well below the random errors. 

The two light curves were then synchronized by finding the 
time shift that minimizes the x 2 of the differences between the 
two data sets, with differences taken between APO data points and 
interpolated McDonald data points. The merged light curve was 
phase-folded according to the period found from the long-term K2 
light curve. The zero phase (defined as the primary mid-eclipse) 
was determined using a parabolic fit around the middle of the pri¬ 
mary eclipse. Fig. 2 shows the combined, phase-folded light curve 
from the fast photometry. We have verified that the light-curve ar¬ 
eas within each eclipse region are consistent between the K2 and 
the ground-based light curves. 

Next, we model the folded light curve during the two eclipses. 
Given the short orbital period of the system, tidal circularisation is 
likely. We will therefore assume a circular orbit, but will test this 
assumption by allowing for a phase shift that differs by 6 from 0.5 
between the primary and secondary eclipses, due to either eccen¬ 
tricity of the orbit or to a Rpmer delay (Kaplan 2010). The eclipse 
duration is related to the system parameters by 


A t 


eclipse : 




- b 2 


( 1 ) 


where P is the orbital period, R\ and R 2 are the primary and sec¬ 
ondary radii, respectively, and a is the semi-major axis. The impact 
parameter is defined SLsb = (a/R\) cos i, where i is the orbital plane’s 
inclination to the line of sight. Since the duration of the eclipse is 
very short compared to the period, a constant relative velocity be¬ 
tween the two components during the eclipses can be assumed. 

To account for limb-darkening we use the limb-darkening law 
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of Claret (2000), 

( 2 ) 

where I is the specific intensity, /jl is the cosine of the angle between 
the line of sight and the direction of the emergent flux from the 
center of the star, and c 1-4 are the limb-darkening coefficients. The 
limb-darkening coefficients for the Johnson-Kron-Cousins UBVRI, 
and for Large Synoptic Survey Telescope (LSST) ugrizy filters have 
been calculated by Gianninas et al. (2013) for various WD effective 
temperatures and surface gravity values. From studying the BG40 
transmission curve and the total instrument throughput of Agile and 
ProEM Camera, we find that a linear combination of the LSST ugr 
filters, 

BG40Sys = 0.1 \u + 0.48g + 0.41r, (3) 

provides a reasonable approximation for the wavelength response 
of the BG40 filter combined with that of the detector and the instru¬ 
ment. For every pair of assumed values of T Q ff and logg, the limb- 
darkening coefficients of each LSST filter were interpolated from 
the table of Gianninas et al. (2013) using bilinear interpolation. In 
our analysis, each WD is modelled as a disc divided into 200 con¬ 
centric annuli. Every annulus is assigned a flux, corresponding to 
the limb-darkening formula on equation (2). The fluxes calculated 
for the different LSST filters are combined using equation (3). A 
model eclipse light curve is then produced by means of numerical 
integration over the uneclipsed regions of each WD. 

The light curve was modelled using a Matlab implementa¬ 
tion 2 of the Goodman and Weare Affine Invariant Markov chain 
Monte Carlo (MCMC) Ensemble sampler (Goodman & Weare 
2010; Foreman-Mackey et al. 2013). We used an MCMC run with 
500 ‘walkers’ and 100,000 steps. The MCMC routine minimizes 
the x 2 value of the fit, over a 5-dimensional parameter space, de¬ 
fined by the ratio of WD radii, R 2 /R 1 , the ratio of semi-major axis 
to primary radius, a/R\, the impact parameter b, the ratio of the 
central surface brightness of the two WDs, / 2 // 1 , and the secondary 
eclipse phase deviation, 6 . The orbital period is taken from the K2 
light-curve analysis. To check for convergence of the MCMC pro¬ 
cess, we verified that the distribution of Ax 2 from the last 75% of 
the steps is well approximated by a x 2 distribution. The medians 
and +1 cr intervals of the parameter distributions in the last 75,500 
steps of the MCMC routine (which are used to estimate the errors 
in the parameters), are listed in the second column of Table 1. The 
third column lists the best-fitting model. 

The limb-darkening coefficients (four for each WD) were kept 
fixed during each MCMC run. The coefficients ( u : 0.92, -0.93,0.92, 
-0.34; g : 0.90, -0.85, 0.72, -0.25; r: 0.85, -0.96, 0.88, -0.32) were 
calculated using equation (3), with the effective temperature and 
logg values found for the primary WD in a spectral analysis ( T\ ~ 
25000K, loggi ~ 7.3, see section 3.3 below), assigned to it. For 
the secondary, T 2 and logg 2 values {J 2 ~ 14000K, logg 2 ~ 7.1, 
derived in Section 3.5 below) were used to evaluate the secondary 
WD limb-darkening coefficients ( u : 0.77, -0.74,0.96, -0.38; g: 0.66, 
-0.13, 0.28, -0.11; r: 0.61, -0.14, 0.07, -0.03). Varying T e ff and logg 
had little effect on the results. 

Table 1 and Figs. 2 and 3 show the results of our modelling 
process. An acceptable best fit is found, with a reduced^ = 1.4. 


2 The code used here is a modified version of Aslak Grinsted’s 
gwmcmc. The unmodified version can be found at https: //github. com/ 
grinsted/gwmcmc 


Table 1 . Eclipse data model results 


Parameter 


MCMC results 

Best-fitting model 
(see Fig. 2) 

R 2 /R 1 


1 14+0.24 

1,1 -0.19 

1.04 

a/Ri 


43 7 +5 - 5 
^-4.1 

41.6 

b 


0 56 +0 - 13 

-0.15 

0.49 

fi/fi 


0.310 + 0.015 

0.310 

5 


0.00008 + 0.00017 

0.00010 

x 2 ■ 

A min 



4414.7 

Degrees of freedom 
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Figure 2. Combined, phase-folded and normalized APO and McDonald 
light curve. The red dashed line represent the best-fitting model (see Ta¬ 
ble 1). Top: full-period light curve; bottom left: zoom on the primary 
eclipse; bottom right: zoom on the secondary eclipse. 


The phase deviation 6 is consistent with zero. There is thus no in¬ 
dication of orbital eccentricity. From the full solution of the sys¬ 
tem parameters (see Section 3.5, below), the expected Rpmer de¬ 
lay is At r0 m = a/c sin/ « 2.5 s, or 6 « 0.00029, which is consis¬ 
tent, within the errors, with our result for 6 . The other parame¬ 
ters are determined to 10-20% precision. The secondary WD has 
a radius comparable to that of the primary, but with a preference 
to be slightly larger. The secondary WD’s lower central surface 
brightness means it is substantially cooler than the primary, as al¬ 
ready clear from the ratio of eclipse depths. The ephemeris for mid¬ 
primary eclipse is BJD(TDB) p = 2457127.801118 + 0.000005 + 
(0.099865 + 0.000013) E, where E is the number of periods since 
UT 2015 April 15,07:05:08. 

We have further searched for the effects of relativistic beaming 
in the K2 light curve. The primary has a radial velocity amplitude 
of K\ « 230kms _1 (see Section 3.5), and its expected beaming 
relative semi-amplitude (Loeb & Gaudi 2003; Zucker et al. 2007) 
is A « BK\/c , with a maximum at phase 0.25 (primary approach¬ 
ing us) and minimum at 0.75 (primary receding), where B is a k- 
correction-like factor depending on the spectrum and on the ob- 
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Figure 3. One- and two-dimensional projections histograms of converged MCMC chains, from the modelling of the eclipse light curves. The contours indicate 
the 1,2 and 3cr levels. Note the strong covariance between R 2 /R 1 , a/R\ and b. The orange markings indicate the best-fitting model (see Table 1). 


served bandpass. Integrating in wavelength over a Planck spectrum 
with the primary WD’s temperature T\ (see Section 3.3) times the 
Kepler response function, following Bloemen et al. (2011), we find 
B = 1.6. For the secondary, with its lower temperature (see Sec¬ 
tion 3.5), we find B « 2.1. Furthermore, the secondary has about 
0.31 the surface brightness of the primary, but (see Section 3.5) a 
surface area larger by (R 2 lR\) 2 , and a radial velocity larger by a 
factor M\/M 2 , but with an opposite sign. Accounting also for the 
uncertainties in all these parameters, the total expected beaming 
semi-amplitude is thus A « 0.00055^ ^^, which could be quite 
small because of cancellation of the beaming effects of the primary 
and the secondary. 

To search for a beaming signal, we have modelled the folded 
K2 light curve as a constant with two half-hour-long ‘boxcars’ 
corresponding to the two eclipses, and modulated the model light 
curve with a sine function with the system’s period and a maxi¬ 
mum at phase 0.25. Varying the sine’s amplitude, the best fit to the 
K2 data is with A = 0.00049 + 0.00019, a marginal 2.5cr detection of 
beaming in the K2 data, within the expected error range (see Fig. 1). 
The eclipses around phases 0 and 0.5, when the beaming effect is 
zero, are brief, and therefore inclusion of the beaming effect does 
not influence the parameters derived above from the ground-based 
eclipse data analysis. We will ignore gravitational lensing effects, 
which are pronounced only during the eclipses, at levels of < 10 -4 , 
well below the precision of the ground-based data. 


3.3 Atmospheric model fit 

To further constrain the system, the spectra were fitted using WD 
atmospheric models (Bergeron et al. 1994; Tremblay & Bergeron 
2009; Gianninas et al. 2011). By experimenting with the models, 
we concluded that there is no evidence for dilution of the B aimer 
lines of the primary by a featureless continuum from the cooler 
secondary, but rather that the secondary strengthens the lines in the 
total spectrum, indicating it is also a DA type. The ultraviolet-bright 
spectral energy distribution (SED) seen in this object’s photometry 
from SDSS and from GALEX (see Fig. 4), indicates a hot primary 
WD. We therefore calculated the expected SED of a hot primary 
DA combined with a cooler DA model spectrum, scaled to have 
0.18 to 0.32 of the total g-band flux in the spectra. This range corre¬ 
sponds to L 2 I (L\ + L 2 ), the fractional contribution of the secondary 
to the luminosity from both WDs, considering the values of / 2//1 
and R 2 /R 1 from the eclipse light curve fit, as well as the different 
limb-darkening profiles in the 4000 - 5000A spectral region of the 
hotter and the cooler WDs. The primary effective temperature, T 1 , 
was searched over a range of temperatures. For each value of T 1 the 
secondary effective temperature, T 2 , was calculated using the mea¬ 
sured central surface brightness ratio and the system response (see 
Sec 3.5). By comparing the models to the observed SED, we were 
able to constrain T\ to ~ 25,500 K, with T 2 ~ 14,300 K. The some¬ 
what poor fit of the GALEX FUV and NUV values (see Fig. 4) 
might be explained by the uncertainty of the Galactic extinction 
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Figure 4. Observed spectral energy distribution for SDSS J1152+0248 
compared to model spectra. Dots, from left to right, are GALEX far-UV, 
GALEX near-UV, SDSS u\ g\ r\ V and z’ fluxes, and UKIDSS Y, J and H 
fluxes. All fluxes have been corrected for Galactic extinction. Models are 
a ~ 25,000 K primary DA WD (dashed line), a ~ 14,000 K secondary DA 
WD (dotted line), and their sum (solid line). 



Figure 5. Atmospheric model fit to the MMT spectrum (a sum of all 
exposures), consisting of a primary DA WD model with T\ = 25,500 K 
and loggi ~ 7.4, contributing ~ 0.8 of the luminosity in this band, plus a 
72 = 14,300 K, logg 2 ~ 7 DA WD model that contributes ~ 0.2. 


values in these bands as estimated by Yuan et al. 2013. Using dif¬ 
ferent extinction values (Seibert et al. 2005; Wyder et al. 2007) did 
not improve the results. Atmospheric model fits to the Baimer lines 
in the WHT and MMT spectra yield similar temperatures, as well 
as providing estimates of the surface gravity of the primary, logg j. 
The fitting results for the primary WD listed in Table 2, are based on 
the sum of all 13 MMT spectra, after correcting all epochs to zero 
velocity (see Fig. 5). The resulting values of T\ and loggi (as well 
as the secondary WD’s T 2 and logg2 from our subsequent analysis, 
see Sec 3.5) were then used to recalculate the limb-darkening co¬ 
efficients (see Sec 3.2), and rerun the MCMC fit. The new results, 
which are those presented in Table 1, were within 1 cr of the results 
of the initial run. 


3.4 Radial-velocity curves 

The kinematic data for SDSS J1152+0248 were modelled as fol¬ 
lows. We measured radial velocities using the cross-correlation 


Table 2. Primary WD model atmosphere fit parameters. 


Parameter 

Value 

Primary spectral type 

DA 

T\ (K) 

25500+1000 

loggi 

7.42 + 0.25 



Heliocentric Julian Date - 2450000 [days] 



Figure 6. Radial velocity data and model fit. Top panel: radial velocity vs. 
time; bottom panel: phase-folded radial velocity curve. The dashed line is 
the best-fitting model. 


package RVSA0 (Kurtz & Mink 1999). We used Hy and higher order 
Baimer lines from the MMT spectra for our velocity measurements, 
with a median error per epoch of lOkms -1 . The measured radial 
velocity (RV) curve was then fit with a sine function, assuming the 
measured K2 photometric period, and with the measured amplitude 
Af meas > and the mean velocity y as free parameters: 

Vr ((p)= 7 + £meas sin (Imp ), (4) 

where p is the orbital phase. The best-fitting model gives K meas = 
109+ 5 km s -1 and y = 55 + 5 km s -1 (see Fig. 6). 


3.5 Combined analysis of system parameters 

We now combine the results from our previous analyses: the period 
P from the K 2 photometry; the ratios R 2 /R 1 , a/R\, b , / 2//1 and 6 
from the fast-photometry eclipse light-curve analysis; the primary 
effective temperature T\ and surface gravity loggi from the model 
atmosphere fit, and the measured RV amplitude AT meas from the 
cross-correlation analysis, to estimate the properties of both WDs 
and the orbital separation (see Table 4). 

In order to derive the secondary WD temperature, T 2 , we com¬ 
puted the ratio between two Planck functions, each integrated over 
the filter plus system response (see section 3.2), for various tem¬ 
peratures and T \, to find the temperature that gives a ratio corre¬ 
sponding to / 2//1 =0.310 + 0.015. This results in a secondary WD 
temperature T 2 = 14350+^ K. 

Next, we tested for the possibility that K meas is being biased 
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low in the cross-correlation analysis by the effect of the (unde¬ 
tected) lines from the secondary WD. To do this, we simulated a 
model spectrum consisting of a hot DA WD combined with a cooler 
DA WD, with luminosity ratios in the range constrained by the light 
curve analysis, and with a range of RV phases for each of the two 
components. Each model spectrum was noised at the 5% level and 
then cross-correlated with a model single hot DA WD as in the real 
observations, for the same orbital phase values as the MMT obser¬ 
vations, to simulate a “measured” RV amplitude. We repeated the 
calculation over a grid of primary and secondary RV amplitudes 
(K\ and K 2 , respectively), for the minimal and maximal values of 
the fractional contribution of the secondary to the luminosity from 
both WDs. The required constraint of a simulated “measured” RV 
amplitude of K me as = 110 + 5kms -1 is shown as the broad blue 
swath in Figure 7. 

To derive the possible mass combinations, allowed by the con¬ 
straints implied by the light curve and the RV curve, we have taken 
a grid of Mi ,2 masses, between 0.2 and 1.4Mo, and calculated the 
corresponding R 12 radii using theoretical WD evolutionary cool¬ 
ing sequences. For a WD of given mass and given core and at¬ 
mosphere compositions at birth, these models calculate the radius 
as a function of the declining temperature. Lacking prior knowl¬ 
edge on the WD core composition and hydrogen envelope thick¬ 
ness, we have considered a range of models: carbon-oxygen cores 
with ‘thin’ (10 -10 mass fraction) or ‘thick’ (10 -4 mass fraction) 
hydrogen envelopes and with a 10 -2 mass fraction helium enve¬ 
lope (Fontaine et al. 2001) 3 , and helium core WDs (Althaus et al. 
2013) 4 . We only considered helium core models with R < 0.08 R©, 
since the densities indicated by the period and the eclipse durations 
imply that both WDs are degenerate. For each combination of core 
compositions of the primary and the secondary, R 12 were estimated 
by interpolating over temperature and mass in the cooling-sequence 
tables, using the derived T\ and T 2 temperatures. (Mi, M 2 ) points 
which gave a radius ratio, R 2 /R 1 , within the limits derived from the 
light curve (see Table 1), were used to delineate the allowed region 
in the (Mi, M 2 ) plane. 

Similarly, to further constrain (Mi, M 2 ), we used Kepler’s law 
to calculate the expected orbital separation over the (Mi, M 2 ) grid, 


a = (A G(Ml+M2 )V (5) 

where P is the orbital period measured from the K2 light curve. 
Again, only grid points with an a/R\ value within the limits derived 
from the light curve (see Table 1), were included in the allowed 
region of the (Mi, M 2 ) plane. 

We then calculated the expected RV amplitudes, over the 
same (Mi, M 2 ) grid, using 


* 1,2 



M2,l 
Mi + M2 


sin i 


( 6 ) 


where the eccentricity, e, was taken as zero, and the inclination, 
i, was taken from the MCMC results (see Table 1), to translate the 
R 2 /R 1 and a/R\ constraints on the (Mi, M 2 ) plane to constraints on 
the {K\ , K 2 ) plane. Combining these constraints with the constraint 
on Kmeas> derived from the cross-correlation simulation, as shown 
in Figure 7, we see that, despite the measured cross-correlation am¬ 
plitude of ~ llOkms -1 , the true RV amplitudes K\ 2 are actually 


3 http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 

4 http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/557/A19 



K\ [km/s] 


Figure 7. Constraints on the RV amplitude (K \, K 2 ) plane. The “allowed” 
area for all possible core compositions is represented by the black con¬ 
tour. The Kmeas = llOkms -1 constraint polygon is derived from the cross¬ 
correlation simulation (see Section 3.5), while the other polygons represent 
the R 2 /R 1 and a/R\ constraints for each core composition combination. 


~ 200 - 250 km s -1 . This can be considered an indirect detection 
of the effect on the lines from the secondary on the total observed 
spectrum. Finally, we transform the allowed region in the (K \, K 2 ) 
plane back to the to the (Mi, M 2 ) plane, using 


Mi, 2 


p(k 1 + k 2 ) 2 

- ^ ^ ■ 7 . - A 2,l 

2 nG sin 3 1 


(7) 


with the results shown in Figure 8. 

Taking into account all possible core compositions, we obtain 
Mi = 0.47 + 0.11M© for the mass of the primary WD, and M 2 = 
0.44 + 0.09 M© for the secondary WD (see Table 3 and Figures 7 
and 8). The combination of a He-core primary WD and a ‘thick’ 
CO-core secondary WD, is the only option which is ruled out by 
the combined constraints. 

The effective temperature and mass were then used to estimate 
each WD’s radius, R 12 , and cooling age, t\£, using the theoretical 
WD evolutionary cooling sequences. The errors on Mi, 2 , Ri, 2 and 
t\ 2 listed in Table 4 include both the statistical errors that prop¬ 
agate from Ti, 2 , and the theoretical uncertainty from the range of 
cooling models included. As seen in Table 3, the model atmosphere 
fitting, combined with the WD evolutionary cooling sequences and 
the constraints derived from the light curve and the RV measure¬ 
ments, indicate both WDs are of low-mass, with Mi ^ 0.47 M© 
and M 2 ~ 0.44M©. Both WDs have radii of ~ 0.02 R©. The cool¬ 
ing age of the primary is « 52Myr, while that of the secondary is 
« 230 Myr. 

The secondary WD radius was recalculated using the radii 
ratio derived from the MCMC, and R\. The surface gravity of 
each WD was then calculated using gi ,2 = GM 12 /R 2 2 : All the de¬ 
rived parameters (R 2 , a , i, T2 , loggi, log# 2 ) were calculated using 
the MCMC chains, while normally distributed chains of the same 
length were created for non-MCMC parameters. The medians and 
lcr intervals of the derived parameters are listed in Table 4. 

The time till merger of the system, due to gravitational wave 
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Table 3. Primary and secondary mass, radius and cooling age from WD evolutionary cooling sequences, constrained by the RV and light curves. References: 
CO core models (Fontaine et al. 2001), He core models (Althaus et al. 2013). 


Primary type 

Secondary type 

Mi (M©) 

m 2 (M©) 

*t (*©) 

R 2 C Rq) 

h (Myr) 

t 2 (Myr) 

He 

He 

0.462 + 0.032 

0.420 + 0.030 

0 0 1 93 +00023 
u.uivjq 0017 

°- 0186 too!3 

6i+ 28 

or 23 

304!^ 

He 

CO, thin 

— 

— 

— 

— 

— 

— 

He 

CO, thick 

0.462 + 0.032 

0.450 + 0.060 

0 0193 +00023 

u.ui ^_ 00017 

0-0217" 

6 i + 28 

Ui -24 

191-52 

CO, thin 

He 

0.416 + 0.046 

0.402 + 0.048 

0 oi94 +00018 
u.uiv M ). 0019 

0-01961°X 

20.3!>« 

279!« 

CO, thin 

CO, thin 

0.421+0.041 

0.389 + 0.031 

0 0192 +00017 
u.uivz_o ooi 8 

05 -H 

© b 

o o 
o o 
+ l 
05 

00 

o 

o 

20 - 6 -40 7 

151+ 62 

1Ji -24 

CO, thin 

CO, thick 

0.412 + 0.052 

0.442 + 0.088 

0 0195 +00019 

u.uivj_o 0020 

0 0220 + ° 0044 
u.uzzu _ 0 0047 

20 - 0 -i°o 3 

1851*“ 

CO, thick 

He 

0.465 + 0.055 

0.414 + 0.036 

0 0205 +00027 

U,U U -0.0027 

0 0189 +0 0018 
u.uiov _ 0 ooi5 

91 c+ 20.2 

zo.J _ 69 

292-57 

CO, thick 

CO, thin 

0.464 + 0.026 

0.412 + 0.018 

0 0203 +00024 

0020 

0 oi79 +0 - 0032 

u.ui /v_o ooi7 

23.6-J™ 

162!g 

CO, thick 

CO, thick 

0.498 + 0.088 

0.444 + 0.066 

0 0195 +00034 
u.uivJq 0033 

0 0219 +0 - 0039 
u.uzi v_ 0 0043 

2 8-n 

186+!® 

Full range 


0.47 + 0.11 

0.442 + 0.088 

0.0197 + 0.0035 

0.0213 + 0.0051 

52 + 36 

230+100 



Figure 8. Constraints on the mass ( M \, M 2 ) plane. Each polygon represent 
a different primary+secondary core composition combination. Note that the 
He + COthin option is ruled out by the constraints. 

losses, calculated using 

5 c 5 a 4 

fmerger “ 256 G3 +M 2 )’ 

is 460+280 M >’ r - 

A distance to the system, 464+®® pc, was estimated using the 
SDSS g-filter magnitude (g = 18.35 mag) and a 0.08 mag Galactic 
extinction correction, compared to the summed g-band luminosities 
obtained from the derived temperatures and radii. 


4 DISCUSSION 

Comparing SDSS J1152+0248 to the five previously known eclips¬ 
ing double WDs (see Table 5 and Fig. 9), we see that all of the 
systems consist of a primary DA (specifically a DAZ in the case of 


Table 4. System parameters 


Parameter 

Value 

Source 

PC h) 

2.39677 + 0.00031 

K2 light curve 

X meas (kms -1 ) 

109 + 5 

RV 

K\ (kms -1 ) 

227 + 26 

Constrained simulated RV 

K 2 (kms -1 ) 

226 + 31 

Constrained simulated RV 

Mi (M©) 

0.47 + 0.11 

Constrained models 

M 2 (M©) 

0.442 + 0.088 

Constrained models 

R X (Rq) 

0.0197 + 0.0035 

Evolutionary sequences 

R 2 (Rq) 

n 09972+0.0064 
U.UZZj _ 0 0050 

Derived 

a (Rq) 

0-86t!? 

Derived 

e 

~0 

Derived 

i(°) 

89 275 +0 - 140 
ov.z / j_o 08 i 

Derived 

T\ (K) 

25500+1000 

SED 

t 2 ( K) 

14350^g 

Derived 

loggi 

7.52 + 0.19 

Derived 

l0gg2 

7.38 + 0.24 

Derived 

h (Myr) 

52 + 36 

Evolutionary sequences 

t 2 (Myr) 

230+100 

Evolutionary sequences 

^merger (Myr) 

460^80 

Derived 

D (pc) 

464!*° 

Derived 

T 0 (BJD(TDB)) 

2457127.801118 + 0.000005 

+ (0.099865 ± 0.000013) E 


GALEX J171708.5+675712, see Vennes et al. 2011; Hermes et al. 
2014) and a secondary WD in a close orbit (with periods ranging 
from ~ 12 min to ~ 6h). In three of the systems, the primary is an 
Extremely Low Mass (ELM) WD (with a mass < 0.25 M©), ac¬ 
companied by a massive WD (> 0.7 M©). In the three remaining 
systems both WD masses are low and comparable, as is the case 
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Table 5. Known eclipsing double WD parameters. Sources: NLTT 11748 (Kaplan et al. 2014), CSS 41177 (Bours et al. 2014, 2015), GALEX J1717+6757 
(GALEX J171708.5+675712, Vennes et al. 2011; Hermes et al. 2014), SDSS J0651+2844 (SDSS J065133.338+284423.37, Hermes et al. 2012), SDSS J0751- 
0141 (SDSS J075141.18-014120.9, Kilic et al. 2014), SDSS J1152+0248 this work. ^Estimated from non-detection of a secondary eclipse in (Kilic et al. 2014). 


Parameter 

NLTT 11748 

CSS 41177 

GALEX J1717+6757 

SDSS J0651+2844 

SDSS J0751-0141 

SDSS J1152+0248 

Type 

DA+D? 

DA+DA 

DAZ+D? 

DA+D? 

DA+D? 

DA+D? 

Mi (Mq) 

0.15 + 0.02 

0.378 + 0.023 

0.185 + 0.010 

0.26 + 0.04 

0.194 + 0.006 

0.47 + 0.11 

M 2 (Mq) 

0.72 + 0.02 

0.316 + 0.011 

>0.86 

0.50 + 0.04 

a Q7+0.06 
' -0.01 

0.442 + 0.088 

Ri (Rq) 

0.0428 + 0.0009 

0.02224 + 0.00041 

0.093 + 0.013 

0.0371+0.0012 

0.155 + 0.020 

0.0197 + 0.0035 

R 2 (Rq) 

0.0109 + 0.0002 

0.02066 + 0.00042 

0.0092 + 0.0026 

0.0142 + 0.0010 

0.0092 + 0.0026 

n 0999+0-0064 
u.uzzo_ 0 0050 

loggi 

6.35 + 0.06 

7.322 + 0.015 

5.67 + 0.05 

6.76 + 0.04 

5.54 + 0.05 

7.52 + 0.19 

log £2 

8.22 + 0.01 

7.305 + 0.011 




7.38 + 0.24 

T\ (K) 

8706+137 

22439 + 59 

14900 + 200 

16530 + 200 

15750 + 240 

25500+1000 

t 2 ( K) 

7594+123 

10876 + 32 

~ 16750 + 3050 

8700 + 500 

< 8900* 

14350+500 

P(. h) 

5.64145164(7) 

2.784370445(36) 

5.907288(72) 

0.212557373(15) 

1.85(99) 

2.39677(31) 

a (Rq) 

1.53 + 0.03 

0.886 + 0.014 




0 .86*!;}? 

i(°) 

89.6 + 0.1 

88.97 + 0.02 

~ 86.8 + 0.6 

84.4 + 2.3 

85.4-+4 ; 2 

89.2753-140 



0 1 - 1 - 1 - 1 - 1 - L - 

0 5 10 15 20 25 

T [10 3 K] 


Figure 9. Masses and temperatures of the components of known 
eclipsing double WDs. The labels mark the primary WD of each 
system: 1. NLTT 11748, 2. CSS 41177, 3. GALEX J1717+6757, 
4. SDSS J0651+2844, 5. SDSS J0751-0141, 6. SDSS J1152+0248 (see Ta¬ 
ble 5 for references). 

in SDSS J1152+0248. CSS 41177, the only double-lined system 
among the six, is almost a twin to SDSS J1152+0248. Both have a 
period of ~ 2h, a primary DA of mass ~ 0.4 M© and temperature 
~ 2 x 10 4 K, and a somewhat lighter secondary with a temperature 
of - 1x10 4 K. 

Badenes & Maoz (2012) have shown that ~ 5% of WDs are in 
short-period ‘double degenerate’ binaries (<a < 0.05 AU ~ IIRq). 
Assuming a separation distribution with equal numbers per loga¬ 
rithmic interval (Badenes & Maoz 2012), the fraction of WDs with 
a < 1%, like SDSS J1152+0248, is ~ 2.5%. Since the probability 
to observe a system such as SDSS J1152+0248 at an inclination 


such that it undergoes eclipses is 

^ 1+^2 /m 
Peclipse — ^ ~ vv 

the probability, per WD, of discovering an eclipsing system is 
~ 0.15%. As mentioned in Section 1, SDSS J1152+0248 was dis¬ 
covered while searching for period variations in the light curves of 
~ 150 WDs in K2 data, among which ~ 125 had signal-to-noise 
ratio sufficient for detecting the eclipse amplitudes seen in this sys¬ 
tem. Thus, we may have been mildly fortunate in discovering one 
such system in the K2 data, but the very-short-period double-WD 
fraction could perhaps be somewhat higher than the few percent 
estimated above. 


5 CONCLUSIONS 

We have discovered that SDSS J1152+0248 is an eclipsing double 
WD, and have used spectroscopy and fast photometry to determine 
the parameters of the system. This short-period system consists of 
a low-mass hot primary DA WD, in a circular orbit with a cooler 
and somewhat lower-mass DA or DC WD with undetected lines, 
but likely also a DA. As in four of the five previously known sys¬ 
tems of this type, the non-detection of spectral features of the sec¬ 
ondary WD prevents a model-independent estimation of the phys¬ 
ical parameters. Nevertheless, these systems permit measuring the 
masses, radii and temperatures of the secondary WD components, 
sharpening our picture of the close double WD binary population. 
Continued K2 observations of WDs in additional fields will likely 
uncover additional eclipsing double WDs, and will thus increase 
the known numbers of these useful systems. 
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